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ABSTRACT 

Confexf. f Pup is the X-ray brightest O-type star of the sky. This object was regularly observed with the RGS instrument aboard 
XMM-Newton for calibration purposes, leading to an unprecedented set of high-quality spectra. 

Aims. We have previously reduced and extracted this data set and combined it into the most detailed high-resolution X-ray spectrum 
of any early-type star so far. Here we present the analysis of this spectrum accounting for the presence of structures in the stellar wind. 
Methods. For this puipose, we use our new modeling tool that allows fitting the entire spectrum with a multi-temperature plasma. 
We illustrate the impact of a proper treatment of the radial dependence of the X-ray opacity of the cool wind on the best-fit radial 
distribution of the temperature of the X-ray plasma. 

Results. The best fit of the RGS spectrum of f Pup is obtained assuming no porosity. Four plasma components at temperatures 
between 0.10 and 0.69 keV are needed to adequately represent the observed spectrum. Whilst the hardest emission is concentrated 
between ~ 3 and 4 R, , the softer emission starts already at 1 .5 R, and extends to the outer regions of the wind. 
Conclusions. The inferred radial distribution of the plasma temperatures agrees rather well with theoretical expectations. The mass- 
loss rate and CNO abundances corresponding to our best-fit model also agree quite well with the results of recent studies of <TPup in 
the UV and optical domain. 
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1. Introduction 

Through their powerful stellar winds, and their violent deaths 
as supernovae, massive stars have a tremendous impact on the 
ecology of their galaxies. The feedback by the stellar winds of 
early-type stars is indeed a major ingredient of the chemical evo- 
lution of the interstellar medium. Moreover, the mass-loss has 
an important impact on the evolution of the star. The quantity 
of material released by the wind as well as its chemical com- 
position depend on the evolutionary stage of the star. A proper 
understanding of the stellar winds of massive stars is therefore 
fundamental for our knowledge of the stellar evolution and the 
impact on the interstellar medium. However, the determination 
of the mass-loss rate of a star is a complicated task. 

Over the last decade, the importance of the structures of the 
stellar winds and their impact on the determination of the mass- 
loss rates became more and more obvious. Observationally, 
these structures manifest themselves through low-amplitude 
variability of emission lines such as Her (Eversberg et al. 1998 ), 
and through difficulties to fit the observed spectra with model at- 
mosphere codes that assume homogeneous winds (Bouret et al. 
2005). These small-scale structures are usually attributed to hy- 
drodynamic instabilities intrinsic to the winds (Feldmeier et al. 
1997). 
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Accounting for the fragmented structure of the winds leads 
to a downward revision of the mass-loss rates by a factor be- 
tween a few and ten. Owing to the different sensitivities of the 
spectra over various wavelength ranges to the structures, multi- 
wavelength analyses are required to determine consistent mass- 
loss rates of massive stars. For instance, Ha emission and free- 
free radio emission are collisional processes which scale with 
the wind-density squared and are thus sensitive to the clumping 
factor regardless of the size of the clumps. In clumped winds 
with a void inter-clump medium, the wind's velocity law misses 
some values along a given sight line. This porosity in the veloc- 
ity field is called 'vorosity' (Owocki l201 11 1. The latter modifies 
the Sobolev length and its impact is mainly visible on the UV 
resonance lines. In a fragmented wind with non-optically thin 
clumps, X-ray lines are sensitive to porosity, which depends on 
clump continuum optical depth, and thus clump size as well as 
the clumping factor. 

The X-ray line profiles can also give information on the hy- 
drodynamical instabilities in the wind. In fact, considering the 
wind embedded shock scenario, the emission of X-ray photons 
arises when two fragments with different velocities collide (Lucy 
& White 1980 ). The observable morphology of an X-ray line de- 
pends on the interplay of emission by the hot plasma and photo- 
electric absorption by the cool wind component. The latter effect 
is obviously also sensitive to the presence or absence of struc- 
tures in the wind. On the other hand, the location of the X-ray 
emission region as well as the temperature of the X-ray emit- 
ting plasma give important information on the hydrodynamical 
properties of the winds. 

£ Pup (= HD 668 1 1 , 04 Ief) has been the target of numerous 
studies over a wide range of wavelengths. This star is often con- 
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sidered as a cornerstone object for the understanding of massive 
stars. The analysis of the far UV - UV - optical spectra by Bouret 
et al. ( 120121 1 confirmed a super-solar total CNO abundance, with 
N being highly overabundant and C being strongly depleted. On 
the other hand, Pauldrach et al. (2012 1 found quite different re- 
sults in their far UV - UV spectral modeling, indicating an over- 
abundance of both nitrogen and carbon and a strong depletion 
of oxygen. Also, they determine a homogeneous mass-loss rate 
(~ 13 x 10 _6 M o yr _1 ) about six times larger than the clumped 
mass-loss rates determined by Bouret et al. (12012) and Najarro 
et al. (1201 11 1. This situation is embarrassing as one would aim 
for consistent parameters obtained from studies using different 
wavebands and different modeling tools. 

£Pup is also a bright X-ray source and has been observed 
both with Chandra and XMM-Newton at high spectral resolu- 
tion. The star was indeed regularly observed with the RGS in- 
strument aboard XMM-Newton for calibration purposes, yield- 
ing a combined data set of unprecedented quality (Naze et al. 
120 1 2 !). Previous studies of £ Pup's high-resolution X-ray spec- 
trum revealed new interesting results (Cassinelli et al. 2001, 
Kahn et al. 120011 Oskinova et al. 12006, Leutenegger et al. 
120061 120071 Coh en et a l. l20T0l i. For instance, th e anal yses by 
Cassinelli et al. (120011 1 and Leutenegger et al. (120061 1 of the 
forbiden, inter — combination, resonance (fir) triplets of He-like 
ions revealed that the inner boundary of the X-ray emission re- 
gion is rather close to the star (at about 1.5 R, and possibly even 
closer for the hottest plasma, Cassinelli et al. 12001;). However, 
these studies were done on individual line profiles considered 
separately. Each specific line was in some sense assumed to be 
formed by a plasma at the temperature of maximum emissiv- 
ity of the line. This approach thus neglects the fact that several 
plasma components (each one characterized by its specific tem- 
perature, emission measure and location inside the wind) may 
contribute significantly to a given line. To the best of our knowl- 
edge, there have been no previous attempts to fit the entire high- 
resolution X-ray spectrum (either RGS or HETG) of an O-type 
star at once by properly combining the contributions from sev- 
eral hot plasma components and modelling the absorption by the 
cool wind, accounting explicitly for the geometry of the prob- 
lem. 

This paper thus explores, for the first time, the possibility 
to fit all the lines of an X-ray spectrum between 6 and 35 A si- 
multaneously and coherently, assuming a discrete distribution of 
plasma components, characterized by their own temperature and 
location in the wind. To do so, we use the code developed by 
Herve et al. (120121 1 which accounts for the intrinsic emissivity of 
the hot plasma as a function of temperature and chemical com- 
position. We further introduce a multi-temperature plasma spec- 
tral fitting procedure. The goal is to determine the properties of 
the X-ray emission region accounting for the most recent results 
from other wavelength domains. These properties can ultimately 
be compared to the predictions of existing and future hydrody- 
namical simulations of the instabilities of stellar winds. 

This paper is organized as follows. In Sect|2] we present a 
short summary of the £ Pup RGS dataset. Then we briefly recall 
the main assumptions in our modelling tool (Herve et al. 120121 1 
concerning the intrinsic emissivity and the absorption formal- 
ism (Sect.O. We also present, in Sect.[4] the different physical 
processes we have included in the code to better reproduce the 
observational spectrum of £Pup. Finally, in Sect.|5] we discuss 
the results that we obtain, especially as far as the composition 
and the location of the different plasma components in the wind 
are concerned. 



2. Observational data 

Naze et al. (120121 1 extracted and reduced 18 observations of 
£ Pup with XMM-Newton which yield more than 700 ks of use- 
ful exposure. For the details of the data reduction, we refer the 
reader to the paper of Naze et al. (120121 Paper I). Using this 
dataset, Naze et al. (2013, Paper II) performed an analysis of 
the variability of the X-ray emission, finding no significant time 
variability of the line intensities and profiles. Therefore, we can 
combine the individual RGS spectra into a single, high-quality, 
fluxed spectrum. The calibration of the fluxed spectrum prop- 
erly accounts for the dependence of the effective area as a func- 
tion of energy, but does not account for the redistribution of 
monochromatic response into the dispersion channels. For in- 
trinsically narrow lines, this might introduce artifacts in spectral 
regions of strong instrumental sensitivity changes. However, in 
the case of £Pup, the intrinsic width of the lines is large and 
therefore such effects should be of very minor impact here. The 
instrumental width of the RGS is accounted for in our fitting 
process, as the synthetic spectra are broadened according to the 
instrumental profile. 

Whilst the Chandra-HETG spectrum of £ Pup (Cassinelli et 
al. l20011 l has a superior spectral resolution and extends to shorter 
wavelengths than the RGS spectrum analysed in the present 
work, it has a considerably lower signal-to-noise ratio than the 
RGS spectrum. We therefore decided to focus our analysis on 
the RGS spectrum only. 

3. X-ray spectral modeling code 

3.1. Basic principles of the code 

In order to analyse high-resolution X-ray spectra of massive 
stars, we have developed our own dedicated code (Herve et al. 
I2012l l. In our model, we assume that the wind hosts several hot 
plasma components. Each plasma component extends between 
an inner and an outer boundary and has a specific kT which is 
constant over that region. The chemical composition is taken to 
be the same for all components. We then discretize the X-ray 
emission region of a given plasma component into a large num- 
ber of cells. We compute the X-ray emissivity of each cell as a 
function of its temperature and the quantity of material inside 
the cell. For this purpose we use version 2.0.1 of the AtomDB 
database (Foster et al. 120121 1 with emissivities computed using 
the APEC code (Smith et al. l200Tlf ! l This code calculates the 
line emission arising from many physical processes such as col- 
lisional excitation, ionization, and dielectronic recombination in 
a collisionally ionized, optically thin isothermal plasma in ther- 
mal equilibrium. The continuum emission is also included under 
the form of the bremsstrahlung emission and the two-photon de- 
cay. 

We assume the plasma to be embedded in the wind and to 
move along with the latter. The velocity of the wind is deter- 
mined by a /Maw i.e. v(r) = v^il - -jf where Voo is the termi- 
nal velocity of the wind and R* is the radius of the star. For each 
cell, we Doppler-shift its contribution into the observer's frame 
of reference accounting for the projection of the wind velocity 
on the line of sight (MacFarlane et al. 1 199 11 >. 

Since the X-ray emitting plasma is assumed to be optically 
thin, an X-ray photon, once emitted, either escapes or is ab- 
sorbed by the cool wind material along the line of sight. Hence, 

1 The interested reader is invited to visit the section 'physics' for a 
precise description of the calculation of the continuum and lines in a 
hot plasma on the AtomDB website at the url http://www.atomdb.org/ 
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the radiative transfer problem reduces to the calculation of the 
optical depth of the fragmented cool wind material, r, along the 
line of sight from the emitting cell to the observer. For this pur- 
pose, we consider that, at the microscopic level, the cross section 
for photoelectric absorption in the X-ray domain is dominated by 
the bound-free transitions from the K and L energy levels. 

If the clumps are optically thick, the effective opacity is de- 
termined by their geometrical properties (Feldmeier et al. 2003 ). 
The effective opacity can be written (Owocki & Cohen 2006): 



K eff 



V_ 

m, 



K 



where / is the scale-length of the clump, m c is its mass, k is the 
mass absorption coefficient and t c the optical depth of a clump. 
This latter can be written: 



/ kM 11 

^^7 = 4^7 



where (p) is the mean density of the stellar wind at position r, / 
the volume filling factor of the clumps, M is the mass-loss rate 
of the star, and v is the velocity of the wind at the radius r which 
is determined by a /Maw (see above). 

Under this assumption, it appears that the effective reduction 
of the optical depth as a result of the fragmentation depends on 
the ratio between the clump scale and filling factor, also called 
the porosity length (Owocki & Cohen 2006). This result can be 
generalized from the optically thick to the optically thin limit 
using a simplified scaled effective opacity (Owocki & Cohen 
2006): 



K eff 
K 



1 



1 + T r 



(2) 



Using this result in a steady state (i.e. non variable) wind, we 
obtain the wind optical depth (Owocki & Cohen 120061) in con- 
ventional (p,z) coordinates: 



r rm 

TAP, Z) t, j _ R ) + ^ ^ h(r/) 



(3) 



where r 



y/p 2 + z 2 , 



_ kM 



- I 



is the characteristic optical 



depth of the wind, and h(r) = j is the porosity length. In the 
following we adopt h(r) — hf x r. 

A slightly different formalism to describe the fragmentation 
of the wind is based on the concept of the fragmentation fre- 
quency (Oskinova et al. 2006). As we have shown in Herve et 
al. (120121 ). both formalisms yield very similar X-ray line profiles 
and the capabilities of the current generation of high-resolution 
X-ray spectrographs aboard XMM-Newton or Chandra do not 
allow us to choose between the two formalisms. Here we adopt 
the porosity length formalism because it leads to an analytic ex- 
pression for the integral in the equation of the optical depth. This 
has the advantage to allow a faster computation. In Herve et al. 
(120121 ). we have shown that the inclusion of a few percent of the 
wind material in a homogeneous inter-clump medium has little 
impact on the emerging X-ray spectrum. In the present work, we 
thus consider that the inter-clump medium is void. 

In either formalism, the characteristic wind optical depth t„, 
and thus the mass absorption coefficient k, play a key role in the 
computation of the absorption of X-rays by the stellar wind. To 
estimate the values of k, we use the non-LTE, line blanketed, 
model atmosphere code CMFGEN (Hillier & Miller [T998j ) to 
compute the ionization structure of the cool wind as a function 



Table 1. Chemical elements and their ionization levels used to 
calculate the mass absorption coefficient k and the UV radiation 
field in CMFGEN. 
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Fig. 1. Radial dependence of the mass absorption coefficient at 
different wavelengths (specified in the upper left corner of each 
panel). The mass absorption coefficients are calculated with the 
elements listed in Tab[T]with abundances found in our study. 



of temperature and mass-loss rate. In these ionization calcula- 
tions, CMFGEN accounts for different physical processes such 
as collisional excitation, photo-ionization by photospheric light 
and X-rays, as well as radiative and dielectronic recombination. 
Once we know the ionization structure, we can then compute k. 

The results are illustrated in Fig.[T] at four different wave- 
lengths. This plot shows that the approximation of a mass ab- 
sorption coefficient that is independent of the position in the 
wind is not a valid one in this case. Indeed, the ionization struc- 
ture of the wind is not constant throughout the entire wind and 
the impact on the value of the mass absorption coefficient is 
huge. The outer regions of the wind are less ionized and con- 
sequently, the cool material produces a stronger absorption in 
the outer regions. 

Compared to our previous version of the code (Herve et al. 
120121 ) where we adopted a value of k that was independent of 
r, we have decided to render the treatment of the absorption of 
the wind more realistic. For this purpose, we approximate the 
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Wavelength in Angstri 



Wavelength in Angstroi 



Fig. 2. Wavelength dependence of the mass absorption coeffi- 
cient at different positions inside the wind (1 .5, 3, 6 and 15 R t as 
specified in the upper left corner of each panel). 



radial dependence of k by two simple relations that hold over 
two ranges of values of r. A linear relation 



K(r) - Ko + y(r - R t ) 

for r < Rum, and a constant 

— K max 



(4) 



(5) 



for r > Rn m . Here kq and K tnax are, respectively, the mass absorp- 
tion coefficient in the innermost part and in the outer part of the 
wind and y is the slope. /?/,,„ is the position in the wind where k 
reaches its plateau value K max (see Fig.[T|). The advantage of this 
parametrization is to provide a more realistic description of the 
variation of k as a function of the radius and, at the same time, 
to preserve an analytical solution for the integral of the optical 
depth equation. Indeed, for a linearly increasing porosity length 
(h(r) — h'x r), the optical depth is now given by: 



ta(p,z) = J 

♦j; 



ch(ko — yR*)R t dz' 



r> (r' -R*) + a(K + y(r' -R*))R»h! r' 



0' - R,) + a (k + y (V - R*)) R* h' 
aK max R*dz' 



V(R/„„) 2 -p 2 r ' (r' ~R*) + a K max R, h' f 



(6) 



with a 



Note that we have assumed in Eq.[6]that z is included in the 

interval [- ^R 2 lim - p 1 , ^jR 2 im - P 2 ]- If me latter condition is not 
fulfilled, the boundaries of these integrals are adjusted, but the 
principle remains the same. 

The first two integrals have analytical solutions making the 
computation easier and the same is true for the last integral 
which is formally identical to the one used in the original pre- 
scription of Owocki & Cohen (2006 see our Eq.[3]). 

Figure|2] illustrates the wavelength dependence of the mass 
absorption coefficient at four different positions in the stellar 
wind. The strongest absorption edge in this figure (at about 26 A) 
is due to N iv. The complex shape of this relation could impact 
some of the frequently used diagnostic tools based on line ratios. 



For instance, the ratio between the flux of N vn Ly a at 24.78 A 
and the flux of the Nvi triplets at 28.8 - 29.5 A, which is used 
in line-by-line analyses as a temperature indicator, could be af- 
fected by the strong N iv edge. 

3.2. The treatment of the fir triplets of He-like ions 

In an optically thin thermal plasma at high temperatures, some 
elements retain only two electrons and are thus present as ions 
with a He-like electron configuration (e.g. O vn and Nvi). Such 
ions produce a triplet of lines consisting of so-called forbidden 
2 3 S\ - 2 'So (f), inter-combination 2 3 P\ t 2 - 2 'So (?) and res- 
onance 2 l Pi - 2 1 So(r) lines (Gabriel & Jordan [T9691 l. In low 
density environments or in the absence of strong UV radiation 
fields, the upper level of the /transition is only depopulated via 
the forbidden line which is then seen as a rather strong line in 
the spectra (Blumenthal et al. 1 1 9721 ). However, in the wind of a 
massive star, the presence of an intense UV radiation field alters 
this situation. Indeed, the UV photons pump the electrons from 
the upper level of the / transition to the upper level of the i line 
(Blumenthal et al. l 19721 Porquet et al. 2001). Therefore the ratio, 
% — I is strongly modified which provides an interesting tool to 
constrain the dilution of the UV radiation field at the position of 
the X-ray emission and hence to determine the location of the 
hot plasma in the wind. 

Since the flux of UV radiation seen by the ions in the stellar 
wind decreases in proportion to the geometrical dilution factor 
w(r) = i(l - [1 - i^fj 1 ] 1 ), the % ratio is also a function of 
radius. The radial modification due to the photoexcitation from 
the metastable upper level of the forbidden line can be expressed 
as follows (Blumenthal et al. l 19721 Leutenegger et al. 2006): 



W) = 



1 



1 + 2Pw(r) + NJN C 



(7) 



where P = £ and d>* = 8tt^ f?& with H v the Eddington flux 

and / is the sum of the oscillator strengths for 2 3 S i to all three 
of 2 3 Pj. N e is the free electron density. (f> c and N c are the critical 
photoexcitation rate and the critical density (which depend on 
atomic parameters and the electron temperature, see Blumenthal 
et al. 119721 ). In the winds of hot, luminous massive stars, the pho- 
toexcitation rate dominates over the collisional excitation term, 
and in the remainder of this work, we have thus neglected the 
NJN C term. 

The treatment of these effects is now included in our code 
along the above prescriptions and using the Eddington fluxes 
computed with our CMFGEN model. In Table|2] we indicate the 
physical information necessary to reproduce this modification of 
intensities of the forbidden and inter-combination line present in 
the spectra of £ Pup. 

4. Multi-temperature plasma fitting procedure 

The full X-ray spectrum of a massive star can usually not be fit- 
ted with a single temperature plasma model. In the present case 
for instance, the most prominent ionization level of iron (Fe xvn) 
seen in the RGS spectrum of f Pup, suggests the presence of a 
hot plasma component with a kT close to 0.50 keV. Other lines, 
such as the Si xm lines around 6.7 A, have their maximum emis- 
sivity at even higher temperatures around 0.90 ke V. On the other 
hand, the fir triplets of He-like carbon and nitrogen are rather 
indicative of a lower temperature plasma (kT around 0.10 keV). 
Fitting the observed spectrum hence requires the combination 
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Table 2. Parameters for the treatment of the forbidden, inter-combination and resonance lines. 



Ion 


X r 




X f 


-t«v(il) 




<f>c 


H v i/(hv) 


H v2 /(hv) 


/ 


Nvi 


28.792 


29.074, 29.076 


29.531 


1907.26 


1896.74 


1.83 x W 1 


1.71 


2.27 


0.1136 


Ovu 


21.603 


21.796, 21.799 


22.095 


1638.28 


1623.61 


7.32 x 10- 


1.65 


1.30 


0.0975 


Neix 


13.447 


13.548, 13.551 


13.697 


1272.81 


1248.28 


7.73 x 10 j 


1.95 


1.87 


0.0700 


Mg xi 


9.1681 


9.2267, 9.2298 


9.3134 


1034.31 


997.46 


4.86 x 10 4 


1.19 


1.10 


0.0647 


Si xiii 


6.6471 


6.6838, 6.6869 


6.7394 


865.14 


814.69 


2.39 x 10 5 


0.64 


0.86 


0.0562 



Notes. /i„,,(,i) and A mii2) correspond to the UV radiation wavelengths which depopulate the inter-combination level of a given element. All wave- 
lengths are given in A. 
<f> c are taken from Blumenthal et al. ( 11972) . 

H v /(hv) are calculated with CMFGEN with our own abundances and mass-loss rate but with T e ff determined in the UV/optical domain. The 
values are given in 10 8 photon cirr 2 s~'Hz~' . 

The oscillator strengths, /, are extracted from the AtomDB website. 



of several models with different temperatures. At this stage, we 
have to stress that the HETG spectrum of ^Pup (Cassinelli et 
al. 120011) as well as its EPIC spectra (Naze et al. 12012.1 ) reveal 
the presence of a weak Si xiv Ly a line at 6. 18 A and of a weak 
S xv triplet around 5 .04 - 5 . 1 A. These ions have their strongest 
emission at temperatures near 1.3 - 1.4keV. However, these lines 
are outside the sensitivity range of the RGS instrument and are 
thus not included in our fit. We will return to this point in Sect. [5] 

Of course, each line is not only emitted at the temperature 
where its emissivity reaches its maximum value, but rather exists 
over some range of temperatures, although with a lower inten- 
sity. Therefore most lines in the observed spectrum result from 
the combination of several contributions associated with plasma 
components at different temperatures (see Fig. 2}. The main dif- 
ficulty in the fitting procedure therefore consists in finding the 
best values of the parameters that describe our synthetic spectra, 
including the filling factor fhotgas (i- e the rat i° of the X-ray emit- 
ting volume and the total volume) for each plasma component. 

We have thus developed a code which combines different 
components at different temperatures in order to reproduce the 
total spectrum of £ Pup by a minimization of the x 2 of the fit. 
We have based our fitting procedure on a multi regression rou- 
tind3 There are quite a number of parameters needed to describe 
a model. These are the general wind and stellar parameters (M, 
Voo, /?, R* and T„), the porosity parameter h' , the abundances of 
the different elements that make up the hot plasma, and for each 
plasma component, its temperature (kT), filling factor ( fhotgas), 
and inner and outer boundary (R,„ and R oiW ). Not all of them are 
treated as free parameters in the present study: only the CNO 
abundances are varied with respect to solar; v m , R t and T„ are 
all taken from Bouret et al. (2012), whilst we adopt /3 = 1 to 
preserve the analytical solution of the optical depth integrals. 
Nevertheless, we are left with 4 x n + 5 parameters where n is the 
number of plasma components. In our case, we find that n — 4 
is required to fit the observed spectrum and we are thus dealing 
with 21 free parameters. In order to obtain the best solution we 
used our synthetic spectra simulator code (Herve et al. 120121) to 
build a database of individual plasma components (hence nine 
free parameters per element of the database). In doing so, each 
parameter is sampled over a given range that we consider rea- 
sonable (e.g. kT varies between 0.08 and 0.80 keV and M varies 
in the range 1 .75 x 10~ 6 to 4.0 x 10" 6 M Q yr _1 ). 



2 We use the IDL routine IMSLJVIultiregress. 
More information are available at 

http://idlastro.gsfc.nasa.gov/idl_html_help/IMSL_Multiregress.html 



The fitting routine then systematically explores this database, 
combining only plasma components that have the same chemical 
composition, h' and M. The synthetic multi-temperature spec- 
trum is built as follows: 

n 

S i = ^ J fhot gas ,jM u + e (8) 
i=i 

where 5, is the observed flux at the wavelength fhotgas, / is the 
classical volume filling factor of the hot gas component j, Mjj is 
the synthetic flux of each X-ray emitting plasma component at 
the wavelength i, and e is the error of regression. 

For each combination of plasma models, the best fit values of 
fhotgasj are determined in the least-squares sense and the result- 
ing xl 1S obtained comparing the model to the fluxed spectrum. 
These numbers are stored for later use, including the identifica- 
tion of the best-fit model based on the value of xl- After sev- 
eral tests, it appears that at least four components are needed to 
achieve a good fit of the fluxed RGS spectrum of ( Pup. 

We have included an additional physical constraint to this 
fitting procedure : we imposed the global X-ray emitting region 
around the star to be continuous. In other words, from the lowest 
inner boundary radius of the emitting regions to the largest outer 
boundary radius, there must not exist any spatial gap where there 
is no hot plasma emission at all. 

It is difficult to make any statement about the unicity of the 
solution. There might indeed exist other combinations of (more 
than four) plasma components with values of the wind param- 
eters (M, h' ,...) outside the range explored here that could fit 
the spectrum. However, within the assumptions made here (four 
discrete plasma components,...), we are confident that our pro- 
cedure has identified the best solutions. 

5. Results and discussion 

To start our analysis, we used the stellar parameters of Bouret 
et al. (120 121 ). Then we explore a large range of values for all the 
free parameters in order to minimize the x 1 of the fit. Our best-fit 
model of the full RGS spectrum of f Pup is illustrated in Fig. [3] 

5.1. Temperature distribution and hot gas filling factor 

The true temperature distribution of the X-ray emitting plasma 
in the wind of ^"Pup is almost certainly a continuum of tem- 
peratures between a minimum and a maximum temperature. 
Nevertheless, a discretization with four different temperatures 
is sufficient to reproduce correctly the observed RGS spectrum 
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Table 3. Best-fit X-ray emitting plasma parameters found in our study. 





model [ 


h'=0.0 


Xv = 15-16 


model 2 


h'=0.02 


Xt = 17.82 


kT 


fhotgas 






fhotgas 


R/« 


Row; 


keV 




R, 


R, 




R, 


R, 


0.10 


0.012 


7.5 


85. 


0.013 


2.8 


95. 


0.20 


0.012 


1.5 


38. 


0.012 


1.6 


40. 


0.40 


0.020 


2.7 


4.0 


0.020 


3.0 


4.2 


0.69 


0.007 


3.1 


4.1 


0.007 


3.2 


4.9 



Notes, x 2 i s the reduced chi-squared with 1222 degrees of freedom. 
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Fig. 3. Our best-fit model (modelj of Table[3] in dashed red line) of the observed RGS spectrum of ( Pup (in solid black line). The 
most prominent emission lines are labelled. We note that there are many more weaker lines that contribute to the spectrum and are 
included in our model. The residuals (in the sense observation minus model) are shown in the panels below the spectrum. 



of £Pup. As pointed out above, these four components might 
not be sufficient to account for the presence of the highly ion- 
ized Si xiv and S xv lines seen in the HETG and EPIC spec- 
tra. However, although the RGS sensitivity range does not cover 
the S xv triplet, our model predicts an integrated line flux of 
0.64 x 1CT 13 ergcrrT 2 s , which is not too far away from the 
value 0.81 x 10~ 13 ergcrrT 2 s _1 reported by Cassinelli et al. 
(2001), especially in view of the relative uncertainties on this 
flux. The latter uncertainties are likely of order 10%, as this 



triplet contains between 80 and 100 counts in the HETG spec- 
trum, as judged from the Chandra transmission grating catalog 
(Huenemoerder et al. 1201 1 PI Therefore, if a hotter plasma com- 
ponent is needed to fully account for the flux of this triplet, its 
emission measure is likely very low. 

In this respect, it is interesting also to compare our results 
with the fits of the EPIC spectra reported in Paper I (Naze et al. 



3 see http://tgcat.mit.edu 
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Fig. 4. Contributions of the four plasma components to the best-fit model of the RGS spectrum of £Pup (modeli, in solid black 
line): kT\ = O.lOkeV (in dotted red line), kTi - 0.20 keV (in dashed green line), kTi = 0.40 keV (in dotted dashed blue line) and 
kT\ = 0.69 keV( in long dashed dotted pink line). 



120121 1. In that paper, we used a four- temperature thermal plasma 
model with an overlying 'slab' of absorbing material. The best 
fit plasma temperatures were found to be kT\ 2j,4 = 0.09, 0.27, 

0. 56 and 2.18keV. Whilst the first three temperatures are very 
similar to - or mean values of - what we find here, the fourth one 
is clearly much higher than found in our analysis. However, the 
emission measure of the 2. 18 keV component in the fits of Paper 

1, is about a factor 100 below that of the 0.56 keV component. 
Therefore, this very hot plasma should indeed have a very lim- 
ited impact on the RGS spectrum analysed here. Nevertheless, 
this remark implies that there exists very hot material in the 
wind of £Pup. It would thus be extremely interesting to collect 
a HETG spectrum of quality comparable to the RGS spectrum, 
and to repeat the kind of analysis presented here, to constrain the 
location of this very hot plasma. 

In previous analyses, the line intensity ratios of hydrogen- 
like and helium-like ions was used as a diagnostic of the X-ray 
source temperature (e.g. Waldron & Cassinelli 2007i. Such ra- 
tios are expected to provide some average temperature between 
the various plasma components that contribute to the line for- 
mation. However, Fig. ^indicates that the relative weights of the 
different plasma components in this average depend upon the el- 
ement that is considered. For instance, in the case of oxygen, the 
O vm Lya line is dominated by the 0.20 keV plasma, whilst the 
Ovn triplet includes contributions from the 0.10 and 0.20 keV 
plasma components. Waldron & Cassinelli (2007) derived kT = 
0.24 keV from the oxygen line ratio, which is consistent with the 
dominant role of the 0.20 ke V plasma in this case. If we consider 
neon instead, we see that the Ne x and Ne ix lines are mostly 
formed in the 0.40 and 0.20 keV components, with minor contri- 
butions from the 0.69 keV component. For these ions, Waldron 
& Cassinelli (2007) obtained kT = 0.33 keV from the line ratio. 



This is consistent with the fact that, in this specific case, the two 
dominant components have their relative importance reverted in 
the lines of the H-like and He-like ions. In summary, whilst we 
confirm that line ratios yield some average estimate of the tem- 
perature, it is impossible to know, from a line-by-line analysis, 
what are the weights of the different plasma components in this 
average. 

The two plasma components with the lowest temperatures 
(kT of 0.10 keV and 0.20 keV) extend over the widest range in 
radius and dominate the emission (model i in Tableland Fig.|4j». 
These emitting regions extend from a region relatively close to 
the surface of the star (7.5 and 1.5 R« for kT of 0.10 keV and 0.20 
keV) to very far out in the wind (85 and 38 R«, respectively). The 
two hotter plasma components (kT of 0.40 keV and 0.69 keV) 
start at 2.7 and 3.1 R, and extend over a smaller range of radii, 
out to 4.0 and 4. 1 R„, respectively. 

With these results, it appears that the region between 3.0 and 
4.0 R* contains X-ray emitting plasma with 3 different temper- 
atures. A priori, in an outwards accelerating wind, one expects 
the velocity jumps due to instabilities to have a lower amplitude 
in the inner parts of the wind. Therefore, one expects that the en- 
ergy available in the shocks is only sufficient to heat the shocked 
matter to rather low temperatures. Farther out in the wind, the 
difference of velocity between two consecutive clumps of matter 
increases and the energy available in the shocks increases ac- 
cordingly. In the outermost parts of the wind, the density and 
hence the emission measures become low and the same applies 
to the velocity jumps since the efficiency of radiative accelera- 
tion progressively drops to zero. 

From this simple reasoning, one would thus expect the 
hottest plasma components to arise from intermediate radii in 
the wind, whilst the softer emission would be produced over a 
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wider ranger of radii starting closer to the stellar surface and ex- 
tending farther away from it. 

With the caveat that we currently ignore where the 2.18keV 
plasma, responsible for the high-energy extension of the EPIC 
and HETG spectra, is located, our results are in qualitative agree- 
ment with these expectations as well as with the 1-D line driven 
instability simulation of Feldmeier et al. ( 119971 1. These authors 
show indeed that shocks can arise very close to the stellar sur- 
face (see Fig. [5] A. Feldmeier, private communication). In our 
best fit model, the temperature of the hottest plasma component 
that is present at a given radial position in the wind, increases 
outwards from the innermost boundary of the emission region 
(1.5 R* where kT reaches 0.20 keV) to a value of kT = 0.69 keV 
in a region between 3 and 4R«. Further out, the temperature of 
the hottest plasma component decreases again. The maximum 
shock temperature hence occurs in the 3 - 4 R* range where the 
velocity jumps predicted by the simulations reach indeed their 
highest values. At positions closer to the stellar surface or far- 
ther out in the wind, the simulations predict velocity differences 
between two consecutive clumps that are smaller, in agreement 
with the lower temperatures found in these regions. 

The low value of the filling factor of the kT = 0.69 keV com- 
ponent indicates that only a fraction of the shocks in the region 
between 3 and 4R„ actually heat the plasma to high tempera- 
tures. To ensure that most of the kinetic energy of a fast moving 
clump which hits a slower clump, moving ahead of it, is indeed 
converted into heat, the collision must occur along the radial di- 
rection. From a statistical point of view, in a 3-D wind, the prob- 
ability that the velocity vectors of two colliding clumps are not 
collinear is most probably larger than the probability that they 
are collinear. Thus only a fraction of the collisions will produce 
the required heating to reach the hottest temperature. 

The overlap between the highest temperature and the lower 
temperature is somewhat inherent to our methodology. Indeed, 
whilst we worked with four different temperatures for the X-ray 
emitting plasma, we consider that the emission region of each 
plasma component is continuous (i.e. there are no gaps) between 
its inner and outer boundary. For instance, our results indicate 
the presence of a &r=0.10keV component relatively close to the 
star to fit the intensity of the lines at wavelengths above 20 A, 
but also in the outer parts of the wind to fit the broad wings of 
the line profiles. In our code we can only assume a continuous 
shell extending all the way from 7.5 to 85 R*. The introduction 
of a discontinuous region would considerably increase the num- 
ber of free parameters in the regression routine. However, the 
scenario of an extended, continuous shell for the coolest plasma 
is supported by two considerations. First, as stated above, non- 
collinear collisions produce cooler plasma, throughout the whole 
wind. Second, in a wind embedded plasma, the hottest plasma 
cools down progressively as it moves outwards with the wind. 
For instance, adopting the parametrization of the cooling func- 
tion used by Feldmeier et al. (1997), we estimate a typical cool- 
ing time of ~ 180 s for the 0.40 keV plasma at the inner bound- 
ary of that component. In the same way, we estimate a cooling 
time of ~ 600 s for the 0.69 ke V component. On the other hand, 
the flow times required for the hot gas to move from the inner 
boundary of its production site (near 2.7 and 3.1 R* for the 0.40 
and 0.69 keV plasma respectively) to a location at 5 R* are es- 
timated at about 16500 s and 13500 s, respectively for the 0.40 
and 0.69 keV gas. Hence the X-ray photons emitted at energies 
of £r=0.10keV and 0.20 keV beyond about 5 R„ exist through 
a combination of direct emission from weaker shocks, and the 
cooling of the hot plasma initially heated to 0.40 and 0.69 keV. 




R in R* 



Fig. 5. Snapshot of the radial velocity in the stellar wind as ob- 
tained in 1-D Line Driven Instabilities calculations (Feldmeier 
2012, private communication). 



5.2. Abundances 

Our best-fit model of the RGS spectrum of £ Pup yields an over- 
abundance of nitrogen along with a depletion of oxygen and 
carbon (see Table|4]i. Our CNO abundances agree qualitatively 
with those of the UV/optical studies of Bouret et al. (2012), but 
are at odds with the result of Pauldrach et al. (20121, at least 
as far as the carbon abundance is concerned. The Lya doublet 
of Cvi at 33.73 A is the only isolated feature of carbon in the 
RGS wavelength domain. Nevertheless, it appears impossible to 
reproduce this line with an over-abundance of carbon as advo- 
cated by Pauldrach et al. ( 120121 1. 

The other elements, that are relevant for the X-ray spectrum, 
are assumed to have solar abundances (Andres & Grevesse, 
1989). We cannot determine the abundances of hydrogen and 
helium directly with the RGS spectrum since there are no lines 
of these elements in the X-ray domain. However, they play a 
key role in the determination of the mass absorption coefficient. 
Therefore, by default, we use the results of Bouret et al. (I2012I I. 

The particular properties of ( Pup can explain the strongly 
non-solar CNO abundances. Indeed, £Pup is a rapidly rotating 
runaway star and its kinematic properties and chemical abun- 
dances could indicate that the star was part of a binary system 
that underwent either an episode of mass transfer and super- 
nova explosion of the companion or a merger event (Vanbeveren 
[20TT1) . 

5.3. Mass-loss rate and fragmentation 

We find a slightly higher mass-loss rate ((3.5 + 0.25) x 
10~ 6 M e yr _1 ) than the UV/optical and IR studie s ((2+ 0.2)xl0' 6 
and 2. 1 x 10 M yr~' respectively, Bouret et al. 120 121 Najarro et 
al. 1201 11 1 but the difference is not important. Moreover, Bouret 
et al. (120121 1 indicate that their mass-loss rate is probably too 
low to be consistent with the line driving of the wind. Our de- 
termination of the mass-loss rate is also in reasonable agreement 
with the value of 2.5 x 10~ 6 M Q yr -1 obtained independently by 
Lamers & Leitherer (11993) from the radio fluxes of f Pup and 
by Oskinova et al. (2007J from a fit of the Ha line and the Pv 
resonance doublet. 

Our results indicate that the wind structure of £ Pup is un- 
likely to be very porous. Indeed, we have tested various as- 
sumptions on the porosity parameter: h' = (no porosity), 
h' = 0.07 (roughly equivalent to the fragmentation frequency 
«o = 1-7 x 10~ 4 s~' derived by Oskinova et al. d2Q06| from a 
line-by-line investigation of the HETG spectrum) and an inter- 
mediate value of h! = 0.02 (model2 in Table[3]l. By far the best fit 
quality is achieved for the model that has no porosity (see Table 
0. A similar conclusion was reached by Cohen et al. (2010 l who 
found that porosity was not required to fit individual lines of the 
HETG spectrum. 
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Table 4. Stellar wind parameters 





T 




yS 


R, 


M 


X(C) 


X(N) 


X(O) 




kK 


km s -1 




Rq 


lO^Msyr" 1 








This work 


40 1 


2300 


1 


17.3 


3.5 


6.00 x 10~ 4 


7.7 x 10~ J 


3.05 x 10~ J 


Bouret et al. 2012 


40 


2300 


0.9 


17.3 


2.00 


2.86 x 10~ 4 


1.05 x 10 -:i 


1.30 x 10~ 3 


Pauldrach et al 2012 


39 


2100 




19 


13.7 


8.2 x 10~ J 


9.13 x 10~ J 


1.14 x 10-3 



Notes. 1 We can not determine the effective temperature of the star with our code. As this parameter is important for the mass absorption coefficient 
calculation, we indicate the value we have adopted in our analysis. The CNO abundances are given as mass fractions. Note that Pauldrach et al. 
( 20 1 2 1 do not assume a /3-law for the wind velocity, but perform a full hydrodynamical computation of the wind. 



From the stellar and wind parameters assumed (R*, v x , (3) or 
derived (M) in this paper, we can estimate what would be the 
threshold value of the porosity parameter h' for clumps to be- 
come optically thick. From Fig.[T| we find that k ~ 170 cm 2 g" 1 
at a wavelength of 20 A and at 4 R„. With the volume filling fac- 
tor of the cool gas (/ = 0.05) derived by Bouret et al. ( 120121 . 
we find that a clump size of / = 0.05 R„, i.e. h' = 0.25 would 
be needed to achieve r e = 1 at 20 A (and at 4R*). A value of 
h' < 0.02, as suggested by our fits, implies that t c = 1 at 20 A 
would only be achieved in the outermost regions of the wind, 
where the porosity has increased significantly (assuming, as we 
have done here, h{r) — h' x r). Therefore, in view of our results, 
it seems that porosity does not play a role in the formation of 
the X-ray spectrum of £Pup. This situation contrasts with the 
results of Oskinova et al. ( 120121 1 for the WN5 Wolf-Rayet star 
WR6. The latter authors concluded that the wind of the WN5 
star had to have a very porous structure to allow the observed 
X-ray emission to escape. 

5.4. The fir triplet formation 

In previous studies (Waldron et al. 120011 Leutenegger et al. 
2006, 120071 Oskinova et al. 120061 . the analysis of the fir triplets 
of the different He-like ions was done line by line in order to 
determine the hot plasma temperature and the dilution of the 
UV radiation field in the X-ray emission zone. For a single 
temperature plasma, this analysis constrains, in principle, the 
position of the X-ray emitting plasma in the wind. 

In our analysis, we clearly show that most spectral lines, in- 
cluding the He-like fir triplets, are the result of the combination 
of emission from several components at different temperatures 
(Figj4ji. This implies that the analysis of the fir triplets, assum- 
ing they are formed by a single temperature plasma and not tak- 
ing into account the abundances of the elements, could be bi- 
ased. Let us illustrate this point by the O vn and N vi triplets. In 
Fig -El one can see that the line profiles and intensities of these 
triplets are reasonably well fitted. However, Fig. 2] shows that 
both triplets are actually formed by contributions of the two low- 
est temperature plasmas which produce quite different line pro- 
files and intensities. In the case of the Nvi triplet, for instance, 
the inter-combination line is stronger than the resonance line for 
the O.lOkeV plasma component (in red in Fig©, whilst the re- 
verse situation applies to the 0.20 keV component. The combi- 
nation of both components fits the observed triplet quite well. 

Figure|6] illustrates the contributions to the observed emis- 
sion of each component of the N vi triplet for a line of sight of 
impact parameter p = 0.5 R». Owing to the pumping by the UV 
photons from the upper level of the /transition, the latter line is 
strongly surpressed at the benefit of the i line. This is especially 



NVI resonance line NVI intercombination line NVI forbidden line 
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Fig. 6. Contribution to the observable emission for the fir com- 
ponents of the N vi He-like triplet. The product of the intrinsic 
emission per solid angle with exp(—T) is shown as a function 
of z (in (p,z) coordinates with p = 0.5 R»). The triplet includes 
contributions from the 0.10 and 0.20 keV plasma components 
starting respectively at r — ^p 2 + z 2 = 7.5 R* and r — 1.5 R*. 



true in the innermost regions of the wind, where the 0.20 keV 
plasma component is the sole contributor to the N vi triplet. As 
a result, Q - - l - and the ratio between the strength of the 
inter-combination line and the resonance line tends towards the 
value of Q predicted by the atomic data. 

We further note the presence of numerous weak lines in the 
AtomDB emissivities. Due to the Doppler broadening of the 
lines and to the limited spectral resolution of the RGS instru- 
ment, these weak lines cannot be seen individually but they con- 
tribute to the total emission. This is particularly important in the 
inter-combination line of N vi. 

As can be seen on Figs.[3]and|7] our model has some of its 
strongest deficiencies in the fir triplets. For instance, the model 
predicts too large a flux in the red wind of the blends of the 
Mgxi and Sixm triplets (Fig. 13, which likely indicates an over- 
estimate of the inter-combination (and possibly forbidden) line 
with respect to the resonance line. Likewise, the resonance lines 
of the N vi and O vn ions are predicted with slightly too strong a 
blue-shift compared to the observations. There are several pos- 
sible reasons for these discrepancies. 

Clearly, the sampling of the emitting plasma into four dis- 
crete temperature components is a simplification that affects the 
results. The actual plasma in the wind of £Pup is more likely to 
have a roughly continuous temperature distribution. 

Another point is the radial dependence of the hot plasma fill- 
ing factor. In our model, we assume that the filling factor of 
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a given plasma component remains constant over the full ex- 
tent of the emission region of this component. This might be an 
over-simplification, as the filling factor could change, especially 
for the 0.10 and 0.20 keV plasma components that span a wide 
range of radii. Although implementing a radial dependence of 
the filling factor into our code is in principle possible, we have 
refrained from doing so, as it would introduce additional free pa- 
rameters, thereby rendering the search for a best fit model even 
more difficult. 

A third avenue to explore could be the velocity-law of the 
hot gas. Throughout this paper, we have assumed that the X-ray 
emitting plasma follows the same /3 = 1 velocity-law as the cool 
gas. However, the hot plasma is highly ionized, thus lacking line 
opacity for the radiative acceleration by the stellar UV radiation. 
Although Coulomb forces help to carry the hot gas along with 
the cool wind, one cannot exclude the possibility that the veloc- 
ity law might differ from the one of the cool wind. Testing the 
impact of the velocity law on our fits is beyond the scope of the 
present paper. 

Finally, more severe discrepancies, of similar type as those 
observed here, motivated Leutenegger et al. (12007) to incorpo- 
rate the effect of resonant scattering into their line-by-line fits 
of the N vi and O vn triplets^ Using the characteristic Sobolev 
optical depth to,* as a free parameter, Leutenegger et al. (2007 1 
were indeed able to improve the quality, of their fits of these 
two triplets, but at the expense of very large, sometimes infinite, 
values of To,*. Such large characteristic Sobolev optical depths 
are problematic though, as they require a filling factor of the hot 
plasma near unity, far away from the values derived here and by 
other studies (Hillier et al. [1993 1. This problem was noted by 
Leutenegger et al. (2007 ), who advocated that such high filling 
factors were only required in small regions where line forma- 
tion takes place. However, such a configuration would then also 
have a strong impact on the formation of other lines, which is 
not seen. This is why we decided not to include resonance scat- 
tering into our model and do not consider this effect as a good 
candidate for solving the observed discrepancies. 

Using the % ratios evaluated on the HETG spectra of £ Pup, 
Cassinelli et al. (2001), estimated a radial range between 1.9 and 
4.0/?, and between 1.7 and 2.25 R* for the formation regions 
of the Mg xi and Si xm triplets respectively. Based on the same 
dataset, Leutenegger et al. (2006 ) found an inner emitting ra- 
dius of the Mg xi and S i xm triplets of (1.43 ± 0.10)/?», whilst 
Oskinova et al. (2006) inferred an inner radius of ^ 1.5 R,. 
At first sight, these results of Cassinelli et al. (2001 1 for the 
Mg xi triplet are in reasonable agreement with the location of 
the 0.40 keV plasma component, which dominates the forma- 
tion of these lines in our model (see Fig.|4j). This good agree- 
ment between the local analysis of Cassinelli et al. (120011 ) is due 
to the narrow radial domain covered by the 0.40 keV compo- 
nent. In contrast, Cassinelli et al. (2001) found that the Ne ix and 
O vii triplets were likely formed below 4 and 10 stellar radii re- 
spectively. At first sight, this result disagrees with the location 
of the 0.20 keV plasma component which extends from 1.5 to 
38/?». However, the actual location of the formation region of 
the observable lines results from the interplay between intrinsic 
emission and absorption, which is rather complex (see Fig. [8). 
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4 Leutenegger et al. ([2007 1 further noted that models without res- 
onant scattering produce an inter-combination line that was not suffi- 
ciently blue-shifted compared to the observations (their Figs. 1 and 2). 
The latter discrepancy is less severe in our model thanks to the com- 
bination of several plasma components with different temperatures and 
locations. 



21 21 5 22 22.5 28 28.5 29 29.5 30 



Wavelength (A) 

Fig. 7. Our best-fit (dashed red line) of four^r triplets present 
in the RGS spectrum of £ Pup (in solid black line). 



For instance, the O vn resonance line has its observable emis- 
sion arising mainly from below 15 R* in not too bad agreement 
with the results of Cassinelli et al. (200T)- 

For the S xv triplet around 5.04 - 5.10 A, Leutenegger et al. 
(2006) obtained an inner radius of only (1.1*2 \)R*, i.e. poten- 
tially much closer to the photosphere than any of the emitting 
regions in our model. However, the HETG spectrum analysed 
by these authors is of rather poor quality at these wavelengths. 
Unfortunately, the RGS instrument does not cover this wave- 
length domain and we cannot verify whether or not the mor- 
phology of this line is well predicted by our model. 



5.5. The Fexvn A 17.096 A line 

The forbidden lines of the He-like ions are not the only spectral 
features arising from metastable levels which are depopulated 
by the strong UV radiation field. Indeed, the transition Fe xvn 
A 17.096 A also arises from a metastable level. Mauche et al. 
(2001 ) have shown that this level can be depopulated by a strong 
extreme-UV continuum radiation between 190 and 410 A. They 
also found that in the case of a stellar surface temperature ~ 
55 kK, the level could be totally depopulated. Since our code 
does not include a full treatment of the atomic transitions of the 
Fe xvn ion and since we have no observational constraints on the 
extreme UV emission of ( Pup, we have included an adjustable 
reduction factor to fit the intensity of this line. Our best-fit results 
correspond to a division of the nominal line intensity by a fac- 
tor 4.5 (see Fig.|9}. Therefore, it seems that the metastable level 
is heavily depopulated, which is not surprising as the effective 
temperature of £Pup is 40 kK. The emission seen in Fig.[9]close 
to 17.1 A actually corresponds to the Fe xvn A 17.051 A line. 

5.6. Additional unsolved problems in the fit 

Beside the discrepancies in the fits of the fir triplets discussed 
above, a few additional regions in the spectra are not well fit- 
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Theoretical plasma emission Observed emission 

in photon/ (s.cm~2.A) " ! in photon/(s.cm " 2.A) 




P/R* = 



P/R* = 




Z/R* 

Fig. 8. Line formation region of the O vn resonance line. The left panels yield the variation of the intrinsic emission per solid angle 
(emissivity multiplied by the density squared) along the z axis of the (p, z) coordinates for different values of the impact parameter 
p. The middle panels provide the optical depth r along the line of sight to the observer (located at z — > °o). The right panels yield 
the contribution to the observable emission, i.e. the product of the intrinsic emission per solid angle with exp(-r). The morphology 
of this distribution results from the combination of the two plasma components at 0.10 and 0.20 keV that extend between 7.5 and 
85 R* for the former and 1 .5 and 38 R* for the latter. 
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Fig. 9. Modification of the theoretical line profile and intensity of 
the Fe xvn A 17.096 A line due to the UV radiative flux. The ob- 
servations are in solid black line, the theoretical prediction with- 
out UV pumping is in dashed green line and our best-fit model 
with an ad hoc reduction of the line due to UV pumping is in 
dotted red line. 



ted. The main one is the region between 17.2 and 17.65 A, 
where the flux is under-estimated in our model (see Fig. [9}- The 
atomic database of the AtomDB model does not contain any 
transitions which could easily explain this strong emission in 
the observation. We have also checked whether there are any 



radiative recombination continua in that spectral range, but we 
could not find any. An instrumental effect can be ruled out as a 
similar bump is seen in the Chandra HETG spectrum of £ Pup 
(Cassinelli et al. l200T1 their Fig. 1) and Waldron et al. (EOOll see 
Fig. 2 in their paper) also find a similar structure in the Chandra 
HETG spectrum of ( Ori. 

The problem could be due to an underestimate of the O vn 
series that affects this part of the spectrum. Although the inten- 
sity of the fir triplet is fitted correctly, an underestimate of the 
strength of the Ovu edge near 20.8 A in our model could ac- 
tually lead to an underestimate of this series. Alternatively, the 
problem could point at a lack of atomic data of some elements 
or some peculiar ionization stages which are complicated to cal- 
culate. The failure to reproduce the flux level in this region has 
also certainly an impact on the determination of the attenuation 
factor of the Fexvn A 17.096 A line as well as on the under- 
estimated continuum beetween 18 and 19 A. We do not think 
that this bump could result from a redistribution of the flux of the 
Fexvn A 17.096 A transition as no such feature is seen in other 
spectra where pumping is very effective (Mauche et al. 1200 11 1. 

We encounter a similar problem in the region between 12.3 
and 13 A and for two more 'lines' at ~ 10.0 and 1 1.0 A. If there 
exists another plasma component, hotter than those included in 
our model, in the wind of £ Pup, some of these features could ac- 
tually be due to numerous weak lines of higher ionization species 
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of iron. However, with our current four-component model, these 
lines are not predicted. 

6. Summary and conclusions 

In this paper, we have presented a multi-temperature embed- 
ded plasma model that fits the RGS spectrum of ^"Pup. Figure[3] 
shows that the intensity of the strongest lines is well repro- 
duced, although there remain some minor differences. As for 
any spectral synthesis code, our method relies on some assump- 
tions that introduce some limitations. For instance, one of these 
limitations is the discretization of the temperature structure of 
the hot plasma in the wind, as we use 'only' four temperatures. 
Moreover, our temperature grid itself is also discretized with a 
step in kT of 0.025 keV. As the line intensities and the fio ratio 
for the treatment of the fir triplets are dependent on the temper- 
ature, the combination of our assumptions in the entire fitting 
procedure can probably explain at least some of the remaining 
small differences between the observed RGS spectrum of £ Pup 
and our best-fit model. Additional aspects to be explored in the 
future are the radial dependence of the hot plasma filling fac- 
tor and the velocity law of the hot gas. We further assumed a 
simplified radial dependence of the mass absorption coefficient 
and used the bridging law of Owocki & Cohen (2006) to calcu- 
late the optical depth from the optically thick to optically thin 
regime. Finally, we stress that, even the best-quality atomic data 
have typical uncertainties of 20%. Some of the differences be- 
tween our model and the observations could thus result from the 
uncertainties on the atomic parameters. 

Nevertheless, our code provides a major step forward in the 
analysis of the X-ray spectra of O-type stars. Indeed, we are now 
in the position to fit the full X-ray spectrum consistently. From 
these fits, we derive the temperature distribution in the wind, 
and constrain the mass-loss rate as well as the abundances of a 
number of elements. 

To the best of our knowledge, our results provide the first 
empirical determination of the temperature distribution of the X- 
ray emitting plasma in the wind of a massive star. In our analysis, 
the shocks producing X-ray emission arise close to the surface of 
the star ( 1 .5 R») and extend over a large zone (up to about 85 R„). 
However the highest temperature (kT- 0.69 keV) is produced in 
a smaller region, at 3 - 4 R*, in agreement with ID hydrodynam- 
ical calculations. 

Another important result of our work is our finding that 
porosity does not improve the quality of the fit of the global 
spectrum. It seems thus that the clumps that make up the wind 
of £Pup must have rather small dimensions. 

Concerning the He-like fir triplets, we have shown that anal- 
yses that neglect the multi-temperature nature of the emitting 
plasma could be biased when attempting to determine the inner 
radius of the emission region of the plasma. This is due to the 
fact that several temperatures contribute to the line emissions. 

We find a depletion of carbon and oxygen and an over- 
abundance of nitrogen. We also determine a mass-loss rate of 
3.5 x 10~ 6 M Q yr _I for £ Pup. These results are in agreement with 
recent determinations by Bouret et al. (2012) and Najarro et al. 
d201 II ). but disagree with the analysis of Pauldrach et al. (1201 21) . 

In future studies, we will apply our code to other O-type stars 
with good quality high-resolution X-ray spectra to generalize the 
results obtained for £ Pup in the present paper. 
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